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We present in-depth measurements of the electronic band structure of the transition metal 
dichalcogenides (TMDs) M0S2 and WS2 using angle-resolved photoemission spectroscopy, with fo¬ 
cus on the energy splittings in their valence bands at the K point of the Brillouin zone. Experimental 
results are interpreted in terms of our parallel first-principles computations. We find that interlayer 
interaction only weakly contributes to the splitting in bulk WS2, resolving previous debates on 
its relative strength. We additionally find that across a range of TMDs, the band gap generally 
decreases with increasing magnitude of the valence band splitting, molecular mass, or ratio of the 
out-of-plane to in-plane lattice constant. Our results provide an important reference for future stud¬ 
ies of electronic properties of M0S2 and WS2, and their applications in spintronics and valleytronics 
devices. 


I. INTRODUCTION 

Transition metal dichalcogenides (TMDs) (MX 2 where 
M = Mo or W and X = S, Se, or Te) are layered semicon¬ 
ducting materials that have recently shown promise for 
application to a wide range of electronic devices, includ¬ 
ing the growing fields of spintronics and valleytronics^. In 
the monolayer limit, TMDs are theorized to have a split 
valence band that is nearly fully spin polarized^ near the 
K and K' points of the hexagonal Brillouin zone. The 
spin polarization is predicted to be reversed between the 
K and K' points, implying existence of an unusual spin- 
valley coupling in the TMDs. The ability to selectively 
populate^“^ these valleys demonstrates their viability for 
use in spintronics/valleytronics devices. 

Despite the importance of the split valence band that 
governs the unique spin and valley physics of TMDs, 
there remain questions regarding the origin of these split¬ 
tings in bulk TMDs. In particular, this splitting is theo¬ 
rized to be entirely a consequence of spin-orbit coupling 
in the monolayer limit, and a combination of spin-orbit 
coupling and interlayer interaction in the bulk limit, but 
there is disagreement^”^ about the relative strength of 
the two mechanisms in the bulk limit. Current electronic 
band structure studies of bulk MoS 2 ^~^^ and WS 2 ^ are 
limited by energy and momentum resolution and lack fo¬ 
cus on the valence band splittings. Experimentally, the 
sizes of the splittings have been characterized, but only 


with limited resolution. 

Here we present high-resolution angle-resolved pho¬ 
toemission spectroscopy (ARPES) data of the electronic 
band structure of bulk M 0 S 2 and WS 2 , and analyze our 
measurements via parallel first-principles computations. 
Our findings resolve previous debates on the role of inter¬ 
layer interaction on the valence band splitting, revealing 
that the splitting in the bulk samples is primarily due 
to spin-orbit coupling and not interlayer interaction. In 
order to gain insight into the nature of this splitting, 
we further investigate its origins, its magnitude, and its 
correlation with the electronic band gap across a range 
of TMDs. Our comparisons for various TMDs reveal 
the interplay between the valence band splitting, band 
gap, molecular mass, and out-of-plane to in-plane lattice 
constant ratio, suggesting how specific TMDs would be 
suited for particular applications. 


II. METHODS 

High-resolution ARPES experiments on M 0 S 2 (WS 2 ) 
were performed at Beamline 10.0.1.1 (4.0.3) of the Ad¬ 
vanced Light Source at a temperature of 150 K using 
50 eV (52-120 eV) photons. The total energy resolution 
was 10 meV (20 meV) for experiments on M 0 S 2 (WS 2 ) 
with an angular resolution (AO) of < 0.2°. The momen¬ 
tum resolution (A/cy) is related to the angular resolution 
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by Afc|| ~ ^2mEk/h‘^ • COS 0 • AO where the kinetic en¬ 
ergy of the ejected electrons (Ek) is dependent on the 
incident photon energy^^. The samples were cleaved in 
situ at 150 K in a vacuum better than 5 x 10“^^ Torr. 

The electronic structures of MX 2 (M = Mo or W 
and X = S, Se, or Te) were computed using the pro¬ 
jector augmented wave method^^’^^ as implemented in 
the package within the generalized gradi¬ 

ent approximation (GGA)^^ scheme. A 15 x 15 x 3 
Monkhorst-Pack k-point mesh was used in the compu¬ 
tations. We used experimental lattice constants and re¬ 
laxed the atomic positions until the residual forces were 
less than 0.001 eV/A. The spin-orbit coupling effects were 
included self-consistently. In order to correct the energy 
band gaps, we also performed calculations with the HSE 
hybrid functional^^’^^. 

Single crystals were synthesized using the Br 2 vapor 
transport technique in a two-zone furnace system. Mo, 
W, and S powders (purity 99.9995%) were mixed in stoi¬ 
chiometric ratios in a sealed quartz tube at 1 x 10“^ Torr 
pressure and annealed to 880 °C for one week to yield pre¬ 
cursor powders. Next, they were loaded into 2 inch diam¬ 
eter tubes with Br 2 liquid, pumped down to 1 x 10“^ Torr, 
and sealed. Typical growth parameters were established 
as heat up for one week from room temperature to 910 °G, 
the low temperature zone was set to 890 °G in one day, 
and the system was kept at this temperature for three 
weeks and controllably cooled down to room tempera¬ 
ture in one week. Synthesized crystals displayed strong 
Raman signals at 384 cm“^ and 407 cm“^ for M 0 S 2 , and 
352 cm“^ and 417 cm“^ for WS 2 . Lastly, monolayers 
exfoliated from M 0 S 2 and WS 2 crystals displayed sharp 
bright luminescence peaks at 1.88 eV and 2.04 eV, re¬ 
spectively. 



FIG. 1. Constant energy maps of the electronic band struc¬ 
ture of bulk (a) M 0 S 2 and (b) WS 2 taken at photon energies 
of 50 eV and 90 eV, respectively. The maps were taken at 
and below the energy corresponding to their respective va¬ 
lence band maximum (VBM) at K. The dashed white line 
indicates the two-dimensional (reduced in the kz dimension 
as shown in the inset) hexagonal Brillouin zone. Evolution of 
the valence band with increase in binding energy about the 
r, K, and K' points can be observed for both samples. An 
energy splitting at the K and K' points can be observed in the 
M 0 S 2 maps, most easily at 170 meV or 260 meV below the 
VBM (white arrows). The splitting for WS 2 is not observed 
170 meV below the VBM, but can be observed 540 meV be¬ 
low the VBM (white arrows). Trigonal warping of the valence 
band about the K and K' points can be observed for both 
samples. 


III. RESULTS AND DISCUSSION 
A. Characterizing the Origins of the Splitting 

Eig. 1(a) and 1(b) show detailed momentum-resolved 
constant-energy maps for bulk M 0 S 2 and WS 2 , respec¬ 
tively. The two-dimensional (reduced in the kz dimension 
as shown in the inset) hexagonal Brillouin zone is over¬ 
laid as white dashed lines. The observed suppression of 
intensity in M 0 S 2 along a specific momentum direction 
near T is due presumably to matrix element effects^^. 
The evolution of the valence band about the T, K, and 
K' points can be clearly observed throughout the range 
of binding energies shown. We offset the two sets of data 
so that the valence band maximum at K is shown at the 
same energy and we compare the two band structures 
as binding energy is increased. The two samples produce 
similar spectra with the most notable difference being the 
appearance of the lower split valence band at K at ap¬ 
proximately 170 meV below the local valence band max¬ 
imum (VBM) for M 0 S 2 • This appears as a small circular 
feature centered at K within a larger concentric feature. 


The analogous splitting for WS 2 becomes apparent in the 
spectra at much higher binding energy, as seen in the map 
at 540 meV below the VBM. The splitting for both sam¬ 
ples is shown by vertical white arrows. Trigonal warping 
of the valence band in the vicinity of the K and K' points, 
similar to that seen in other bulk TMDs^’^^, is observed 
for both samples and presented here for the first time for 
bulk WS 2 . This deviation from a circular feature about 
K and is more pronounced at higher binding energies 
where the feature becomes more triangular. An analo¬ 
gous trigonal warping has been predicted^^’^^ to occur in 
monolayer M 0 S 2 samples as a consequence of the orbital 
structure of these bands and anisotropic (order q^) cor¬ 
rections to the energy but has not been observed so far 
in an ARPES experiment. 

To obtain information about the kz dispersion, i.e. 
the interlayer dispersion, in Eigure 2(a) we show de¬ 
tailed photon energy dependent maps compared with 
kz-resolved theoretical calculations of bulk WS 2 along 
the T-K high symmetry direction. Indeed, at normal 
emission {k\\ =0), kz is related to the kinetic energy of 
the ejected electrons (Ej^) by kz = y/2m(E^^rV%J7^ 
where Vin is the inner potential determined from the 
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FIG. 2. (a) Photon energy dependence from 52 to 120 eV 
and theoretical calculations of the electronic band structure 
of bulk WS 2 along the F-K high symmetry direction. The 
horizontal red (yellow) dashed lines correspond to the calcu¬ 
lated location of bulk WS 2 bands at the F point (A point). 
Observed intensity shifts of the valence band near the F point 
at different binding energies is evidence of strong kz disper¬ 
sion. (b) Energy distribution curves (EDCs) (solid black 
lines) taken along k^; = k^^ = 0 of the respective maps in 
(a) (e.g. along the white dashed line in the 52 eV map). The 
vertical red (yellow) dashed lines correspond to the calculated 
location of bulk WS 2 bands at the F point (A point). The 
peak in intensity (denoted by the tick marks) shifts period¬ 
ically between higher and lower binding energies as a result 
of the strong kz dispersion. Peaks in the 70 eV and 110 eV 
EDCs as well as the top of the feature in the 90 eV EDO 
correspond to the theoretical band maximum at A. 


energy (up to 880 meV) of the bands as kz changes from 
0 to ^. Similar kz dispersions are also found theoretically 
and experimentally more generally in TMDs^^~^^. 

Near K, we observe the band splitting at the top 
of the valence bands for all photon energies, indicating 
that the splitting occurs for all corresponding kz values. 
We measure the magnitude of the splitting to vary be¬ 
tween 414 meV and 441 meV in accordance with our 
corresponding predicted range of 410 meV {kz = ^) to 
466 meV {kz = 0). This relatively small change in magni¬ 
tude of the splitting (27 meV experimentally and 56 meV 
theoretically) across this range of photon energies indi¬ 
cates weak kz dispersion near K and K^ This is in con¬ 
trast to the observation of large variations in the spectra 
of the broad feature at F which indicates strong kz dis¬ 
persion of the bands near F. 

As a result of the layered structure of bulk TMDs, 
the three-dimensional character of the band structure is 
directly related to the strength of the interlayer interac¬ 
tion. The strong kz dispersion at F is a reflection of the 
three-dimensional character of those bands. However, 
the relatively weak kz dispersion at K and K' is a reflec¬ 
tion of the two-dimensional character of those bands and 
shows that interlayer interaction has a weak effect on the 
valence band splitting. Hence, we conclude that the split¬ 
ting is dominated by spin-orbit coupling with interlayer 
interaction playing a minor role. 


measured dispersion^^’^^. Near F, we observe a broad 
feature likely associated with the convolution of the var¬ 
ious bands predicted for different kz values as it extends 
between the two horizontal red dashed lines correspond¬ 
ing to the calculated energies of the top two valence bands 
at F {kz = 0). The top of this broad feature aligns well 
with the calculated valence band maximum 8it T {kz = 0) 
for the 52 eV map and as the photon energy is increased 
further, the top of the feature decreases in energy until 
it aligns with the calculated valence band maximum at 
the A point {kx = ky = 0, kz = ^) indicated by the hor¬ 
izontal yellow dashed line in the 90 eV map. Increasing 
the photon energy further returns the top of the valence 
band feature to its maximum value. 

The observed kz dispersion at F is investigated more 
closely in Fig. 2(b) by comparing energy distribution 
curves (EDCs) taken at F for each of the incident photon 
energies (e.g. along the dashed white line in the 52 eV 
image in (a)). The observation of a large spread in en¬ 
ergy of the main peak in the EDC spectra between the 
calculated energies of the top two valence bands at F 
(vertical red dashed lines) is in agreement with the asso¬ 
ciation of the broad feature near F with the convolution 
of multiple bands predicted for different kz values. The 
peak intensity (tick mark) moves up and down in bind¬ 
ing energy in a periodic manner confirming the strong kz 
dependence of these bands. Overall, our observed kz dis¬ 
persions agree well with other predictions for WS 2 ^ and 
our calculations which predict wide variation in binding 


B. The Effect of Spin-Orbit Coupling on the 
Valence Band Structure 

Figure 3 shows high-resolution energy versus momen¬ 
tum maps of the band structure of bulk M 0 S 2 (panel a) 
and WS 2 (panel b). Band dispersions can be seen as far 
as 8 eV below the Fermi energy; many of these bands 
have not been observed before. The experimental band 
structure shows remarkable agreement with our theoret¬ 
ical predictions (blue dashed lines) throughout this en¬ 
ergy and momentum region. The intensity of the valence 
band maximum near F for M 0 S 2 is reduced due presum¬ 
ably to matrix element effects^^. Additionally, intensity 
of the experimental bands with low binding energy may 
be shifted near F due to their strong kz dependence (as 
shown for WS 2 in Fig. 2). The splitting of the valence 
band at K is fully resolved for both samples (as indicated 
by the white arrows) and is found to be 170 ± 2 meV for 
M 0 S 2 and 425 ± 18 meV for WS 2 (as shown in the upper 
insets). The value for M 0 S 2 agrees well with our exper¬ 
imental constant energy maps presented in Fig. 1(a) as 
the lower split band for M 0 S 2 appears 170 meV below 
the higher split band. The value for WS 2 is consistent 
with our experimental constant energy maps in Fig. 1(b) 
as the lower split band is not visible 170 meV below 
the higher split band, but is clearly visible 540 meV be¬ 
low. These measured magnitudes are smaller by approx¬ 
imately 100 meV than the values of 265 meV for M 0 S 2 
and 570 meV for WS 2 predicted by HSE-based calcula- 
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FIG. 3. High resolution maps of the electronic band structure 
of bulk (a) M 0 S 2 and (b) WS 2 along the F-K high symme¬ 
try direction with theoretical calculations overlaid as dashed 
blue lines. Spin-orbit coupling and interlayer interaction in¬ 
duced splitting in the valence band near the K point is fully 
resolved for both samples as indicated by the white arrows 
and the respective EDCs (upper insets) taken along the up¬ 
per white dashed lines. [Here the WS 2 EDC was integrated in 
momentum between the two white tick marks.] The splitting 
is found to be approximately 170 ± 2 meV for M 0 S 2 and 425 
± 18 meV for WS 2 . An additional splitting in a higher bind¬ 
ing energy WS 2 valence band is observed within the white 
dashed rectangle in panel (b). The lower inset presents an 
EDC taken along the vertical white dashed line correspond¬ 
ing to the energy and momentum location of the maximum 
splitting of 196 ± 22 meV. [Here the EDC was integrated in 
momentum between the two white tick marks and smoothed 
with a Savitzky-Golay filter.] The calculated band structure 
is found to be in excellent agreement with the experimental 
band structure for both samples. The discrepancy in the lo¬ 
cation of the valence band near the F point for bulk M 0 S 2 
can be attributed to its strong kz dispersion. 


tions. However, our experimental values agree very well 
with GGA-based calculations that predict a splitting of 
145-224 meV for M 0 S 2 and 410-466 meV for WS 2 across 
all kz values. The large difference in magnitude between 
the splittings for M 0 S 2 and WS 2 can be explained by the 
larger intrinsic spin-orbit coupling of tungsten compared 
with that of molybdenum. Additionally, our calculations 
predict a substantial splitting of a valence band of WS 2 at 
high binding energy (see white dashed rectangle). The 
maximum splitting predicted is in agreement with the 
experimental value of 196 ± 22 meV extracted from the 
EDC spectra (bottom inset) where two clear peaks can 
be distinguished. Similar to the splitting in the top va¬ 
lence band, we predict this band to be split largely due 
to spin-orbit coupling and only partially by interlayer in- 




Splitting (meV) (eV) 


EIG. 4. (a) Magnitude of valence band splitting at the K 
point for bulk MoS 2 ®’^^ MoSe 2 ^^’^^ MoTe 2 ^^’^^ WS 2 ^ and 
WSe 2 ^^’^^. (b) Magnitude of the splitting vs. band gap (Eg) 
for bulk MoS 2 ^^’^®■^^ MoSe 2 ^®’^®■^^ MoTe 2 ^®’^^ 
and (c, d) Ratio of out-of-plane to in-plane 

lattice constant (y)^^’^^ vs. magnitude of valence band split¬ 
ting (c) and band gap (d). The dashed lines are guides to the 
eye indicating the different trends for compounds containing 
molybdenum (red) and those containing tungsten (blue). The 
trends are further subdivided between experimental (circles) 
and theoretical (triangles) data. Eor all plots, filled (open) 
symbols represent data from this (other) work(s). 


teraction. 


C. Relationship between Band Gap and Spin-Orbit 
Splitting 

Figure 4(a) shows the magnitude of valence band split¬ 
ting for bulk M 0 S 2 , WS 2 , and other bulk TMDs. The 
molybdenum- and tungsten-based TMDs are ordered 
separately from left to right with increasing molecular 
mass. For a given transition metal, the magnitude of the 
splitting increases as molecular mass increases. These 
trends are shown by the dashed lines and can be ex¬ 
plained as a consequence of the larger intrinsic spin-orbit 
coupling of the successively heavier atoms. The same 
general trend is obtained from HSE hybrid functional cal¬ 
culations (triangles), although the splitting appears to be 
approximately 100 meV larger. Figure 4(c) shows the ra¬ 
tio of out-of-plane to in-plane lattice constant (^) versus 
the magnitude of valence band splitting. Similar to the 
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molecular mass trend, we observe that for a given transi¬ 
tion metal, the magnitude of splitting increases with in¬ 
creasing ^ ratio; however, our calculations indicate that 
although this trend holds across a range of TMDs, it does 
not hold for individual TMDs. 

In addition to the valence band splitting, the band 
gap plays a key role when considering TMDs for use 
in spintronics/valleytronics devices. Figure 4(b) shows 
the magnitude of valence band splitting versus band gap 
(Eg) for bulk M 0 S 2 , WS 2 , and other bulk TMDs. For a 
given transition metal, the band gap of bulk TMDs de¬ 
creases with increasing valence band splitting. The same 
trend is observed as a function of molecular mass in that 
the band gap decreases with increasing molecular mass. 
Again, this trend is seen for both experimental and the¬ 
oretical data even though the magnitude of the splitting 
for each TMD is calculated to be approximately 100 meV 
larger than the corresponding experimental values. Fig¬ 
ure 4(d) shows the ratio of out-of-plane to in-plane lat¬ 
tice constant versus band gap. It too follows the same 
trend as molecular mass and valence band splitting in 
that the band gap decreases with increasing ^ ratio. This 
strong relationship between the band gap, the valence 
band splitting, and ^ ratio would help guide the search 
for suitable TMDs targeted at specific applications. 

D. Conclusions 

We have presented in-depth high-resolution ARPES 
studies of the electronic band structure of bulk M 0 S 2 
and WS 2 together with detailed first-principles compu¬ 
tations. We have shown that the valence band splitting in 
bulk WS 2 is primarily due to spin-orbit coupling and not 
interlayer interaction, resolving previous debates on their 
relative role. These results strengthen the connection be¬ 
tween bulk and monolayer TMDs as the splitting in the 


monolayer limit is due entirely to spin-orbit coupling. We 
have shown how the magnitude of this splitting changes 
for various TMDs and how it correlates with molecular 
mass, ^ ratio, and band gap. These results provide im¬ 
portant new information on how to control these param¬ 
eters in TMDs and for the development of TMD-based 
devices in the fields of spintronics and valleytronics. 
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